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Abstract 

This paper investigates instability ranges of a cylindrically sym¬ 
metric collapsing cosmic filamentary structure in Brans-Dicke theory 
of gravity. For this purpose, we use perturbation approach in the 
modified field equations as well as dynamical equations and construct 
a collapse equation. The collapse equation with adiabatic index (T) 
is used to explore the instability ranges of both isotropic as well as 
anisotropic fluid in Newtonian and post-Newtonian approximations. 

It turns out that the instability ranges depend on the dynamical vari¬ 
ables of collapsing filaments. We conclude that the system always 
remains unstable for 0 < T < 1 while T > 1 provides instability only 
for the special case. 
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1 Introduction 


Dark energy and gravitational collapse are the most fascinated and inter¬ 
esting phenomena of cosmology as well as gravitational physics. Number of 
astronomical observations such as Supernova type I, Sloan Digital Sky Sur¬ 
vey, large scale-structure, Wilkinson Microwave Anisotropy Probe, galactic 
cluster emission of X-rays and weak lensing describe accelerated behavior 
of the expanding universe [Tj. It is suggested that a mysterious type of 
energy known as dark energy is responsible for this accelerated expansion 
of the universe. This induces the problem of correct theory of gravity and 
thus numbers of modihed theories of gravity are constructed using modihed 
Einstein-Hilbert actions. Brans-Dicke (BD) theory is one of the most ex¬ 
plored examples among various modihed theories that provides convenient 
evidences of various cosmic problems like inhation, early and late behavior of 
the universe, coincidence problem and cosmic acceleration [2]. This is a gen¬ 
eralized form of general relativity (GR) which is constructed by the coupling 
of scalar held 0 and tensor held R. It contains a constant coupling param¬ 
eter ojbd (tuneable parameter) which can be adjusted according to suitable 
observations. This theory is compatible with Mach’s principle, weak equiva¬ 
lence principle and Dirac’s large number hypothesis [3]. It is also consistent 
with solar system observations and experiments (weak held regimes test) for 
|a;| > 40,000 @]. 

Gravitational collapse is a process in which stable stellar objects turn 
into unstable ones under the ehects of their own gravity. The formation and 
dynamics of large scale structures such as stars, celestial cluster and galaxies 
are investigated through this phenomenon. It is believed that diherent in¬ 
stability ranges for astronomical bodies lead to diherent structure formation 
of collapsing models. Ghandrasekhar [5] was the hrst who explored stability 
ranges of a spherically symmetric isotropic huid in GR. He used equation 
of state involving adiabatic index (T) and concluded that the huid remains 
unstable for T < |. Later on, many researchers BE] investigated dynam¬ 
ical instability of diherent types of huids (anisotropic huid, adiabatic, non 
adiabatic as well as shearing viscous huid) in spherical as well as cylindrical 
conhgurations and found that stability ranges depend on physical properties 
of the respective huid. 

It is believed that the study of collapse phenomenon in modihed theories 
may reveal modihcation hidden in the formation of astronomical structures 
[H]- In 1969, Nutku [9] explored instability ranges of spherically symmetric 
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isotropic fluid in BD theory and concluded that BD fluid remains unstable 
for r > |. Kwon et ah m discussed instability analysis of the Schwarzschild 
black hole in BD gravity. Sharif and Kauser [m investigated stability ranges 
for spherical as well as cylindrical collapsing models in f{R) theory and found 
that instability ranges depend upon characteristics of fluids and dark energy 
components. Sharif and Yousaf [12] studied the effects of electromagnetic 
held on instability ranges for various models of f{R) gravity. Sharif and 
Rani [T3| explored dynamical instability of spherically symmetric huid in 
f{T) theory and concluded that modihed terms control instability ranges. 
In a recent paper [Tl|, we have discussed collapse of spherically symmetric 
anisotropic BD huid through instability analysis and found that 0 < B < 1 
always leads to unstable conhguration while B > 1 provides instability only 
for one particular case. 

The behavior of hlamentary structures has important implications for 
the formation of structure in the universe. Galaxy hlaments are the largest 
known cosmic structures in the universe. The hlamentary structure is always 
present in the interstellar medium and instabilities within these hlaments 
create dense medium (dense core) where stars form [TJ] . N-body simulations 
of formation of large scale structure describe a wide range of hlaments (with a 
cluster of galaxies forming at the intersection of hlaments) [T6|. Filamentary 
structures are associated with cosmic web on the large scales and are used to 
describe tidal tails (thrown oh by merging galaxies) on the small scales [T^ . 
In order to understand fragmentation of hlament structures, cylindrically 
symmetric hlament models are widely studied [T8] . 

In this paper, we investigate dynamical instability of cylindrically sym¬ 
metric hlaments collapsing structure in BD gravity. The paper is organized 
in the following format. The next section discusses BD equations, Darmois 
junction conditions and dynamical equations. In section 3, we use perturba¬ 
tion technique to construct hydrostatic equilibrium (collapse equation) and 
describe instability ranges (at Newtonian and post-Newtonian (pN) limits) 
for isotropic as well as an anisotropic huid distributions. Finally, the last 
section summarizes the results. 
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2 Brans-Dicke Theory and Dynamical Equa¬ 
tions 

The BD theory (with self-interacting potential V (0) ) has the following action 

0 

S = J d^x^[cj)R - - V{<P) + Lm], (1) 

where SttGo = c = 1 and Lm represents matter distribution. Varying Eq.(IT]) 
by Qa/B and (p, we obtain the following BD equations 

+ ( 2 ) 

0(1, = ——— +-1- -2V((t,)]. (3) 

3 -|- 2ubd 3 -|- 2ojbd dcp 

Here Ga/B is the Einstein tensor, is the energy-momentum tensor for 
matter distribution with T™ as its trace and □ represents d’Alembertian 
operator. The energy distribution due to scalar held is given by 

/.</>’''] “ ^^^^9afB- (4) 

Equation (|2]) gives the BD held equations and ([3]) is a wave equation for the 
evolution of scalar held. 

We split 4D geometry into interior and exterior regions by considering 
a timelike 3D hypersurface as an external boundary of the respective 
cylindrical body [7]. The interior region of a collapsing cylindrical hlamentary 
structure is represented by 

r)dt^ — r)dr‘^ — r)d(j)^ — dz‘^, (5) 

where we consider comoving coordinates inside the hypersurface. In order 
to preserve cylindrical symmetry, the coordinates satisfy the following con¬ 
straints 

— oo < t < oo, 0 < r < oo, —oo < z < oo, 0 < (p < 2 tt. 

In stationary or static region, a scalar held becomes constant and all station¬ 
ary black holes in BD gravity are identical with GR solutions [IH]. Therefore, 
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for exterior region to we take line element of static cylindrical black hole 
given by 

+ 2drdy — r‘^{d(j)^ + ( 6 ) 

where M, v and 7 describe the total gravitating mass, retarded time, and ar¬ 
bitrary constant, respectively [20]. The interior region is hlled with anisotropic 
matter distribution represented by 

= (P + Pr)UaUp - PrQap + (Pz “ Tr)5'«S'y3 -1- {p^ - Pr)KaKp, (7) 

where p, Pr, Ptj, and p^ indicate energy density and principal pressure stresses, 
respectively. The four velocity Mq, unit four-vectors Sa and are calculated 
as Ua = AS^, Sa = do and Ka = C6'^ satisfying u°'Ua = 1 , S^Sa = K^Ka = 
— 1, S^^Ua = K°‘Ua = S'^Ka = 0. For the interior region, the BD equations 
are given in Appendix A. 

Junction conditions provide smooth connection between interior and ex¬ 
terior regions over We consider Darmois junction conditions to discuss 
connection between two regions [7] and for this purpose we take C-energy 
(mass function) [ 21 ] given by 

E(t, r) = m{t, r) = -(1 — ^fVj^f). ( 8 ) 
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Here E{t, r) is the gravitational energy per unit specihc length of the cylinder, 
r represents the areal radius, p shows the circumference radius and I indicates 
specihc length. These are given as follows 

r pi, p ^ ^(2)/9^(2)’ 

where ^(i) = ^ and ^( 2 ) = ^ are the respective Killing vectors. For the 
interior spacetime, Eq.(( 8 ]) takes the form 


m(t, r) 




(9) 


where dot and prime show derivatives with respect to t and r, respectively. 
Since in BD gravity, scalar held and metric tensor are indicated as gravi¬ 
tational variables, therefore (f) = 0 s(=) = constant at the hypersurface 
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The continuity of first and second fundamental forms (Darmois conditions) 
yield the following relations 


^(e) / 

r = r^ce) = constant, ni(t,r) — M = -, I = AC, 


Pr —Til Ti 


01 


AB 


V{4>) 

24 , 


( 10 ) 


Dynamical equations obtained from the contracted Bianchi identities describe 
the conservation of total energy of the system given by 




q;/3 


+ 


Up = 0 , 


rpaP 

^ m _|_ 4> 


Xp = 0 , 


( 11 ) 


\ / ;Q( \ / ;a 

where X/3 = —BS^ (unit four-vector) which provides 

( 12 ) 

(13) 

Ki and K 2 are mentioned in Appendix A. 


P_ 

A 


B C 

^ + (p + p0jg + (i> + p*)^ 


P'r 4>'Pr / \ 4 \ C" 

B + ^ BC, 


-l- ATi — 0, 

+ K2 = 0, 


3 Instability Analysis 

Here, we use perturbation approach to construct collapse equation which will 
be used for instability analysis. We assume that initially, the system is in 
static equilibrium (metric as well as material parts have radial dependence 
only) and after that all the dynamical variables along with metric functions 
are perturbed and time dependence appears [7j. The scalar held, scalar po¬ 
tential and metric tensors have the same time dependence, while the density 
and pressure bear the same time dependence as follows 


A{t,r) 

= Ao{r) + eT{t)a{r), 

(14) 

B{t,r) 

= Bo{r) + eT{t)bir), 

(15) 

C{t,r) 

= Co{r) + eT{t)c{r), 

(16) 

(j){r, t) 

= (j)o{r) + €T{t)^{r), 

(17) 
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Pr{t,r) 

= Pro{r) + epr{t,r), 

(18) 

Pcp{t,r) 

= P<po{^) + (^P<p{^,r), 

(19) 

p{t,r) 

= po{r)+ ep{t,r), 

(20) 

U(0) 

= Vo{r) + eT{t)V{r), 

(21) 


where 0 < e <C 1 and the static distribution is expressed by zero subscript. 
For static and perturbed conhgurations of the held as well as dynamical equa¬ 
tions, we take Cq = r. The static conhguration of BD formalism, perturbed 
form of BD equations and junction condition flTUD are given in Appendix A. 
The perturbed distribution of hrst Bianchi identity gives 


P = 


{PO+Pro)b (po+P0o)c <hpo , A t7 
H-h —-h AocpoKi 


Bo r 

The perturbed form of Eq. (lT3ll provides 

P'r + {P + Pr)~^ + {Pr “ = 0 , 

Ao r (poBo 


T. 


( 22 ) 


(23) 


where Ki and K 2 are given in appendix A. Equation fl 5 ^ along with junction 
conditions provides 


T(f) = Cl exp(7s(e)f) C2 exp(As(e)f), 


(24) 


where 7 s(e) = +i/f, = -y/f with 


s(=) *F 


2 c 

rAg’ 





00 


1 

Bor 


and Cl, C 2 indicate arbitrary constants. Equation fl2T)) shows static and 
non-static distributions leading to stable as well as unstable phases of grav¬ 
itating system. For instability analysis, we assume that when the instabil¬ 
ity phase begins, the system was in complete hydrostatic equilibrium, i.e., 
(t = —00, r(—00) = 0 ). Using this assumption in Eq. (l2T|) . we have C 2 = 0 
whereas Ci = — 1 is chosen arbitrarily. The corresponding result is described 
by 


T{t) = -exp( 7 s(e)f). 

For a real instability regime, we assume only positive values of A 


( 25 ) 
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For the investigation of instability ranges, we use an equation of state 
involving adiabatic index F [22] given by 


Pj 


F- 


Pjo 


po + PjO 


p- 


( 26 ) 


The adiabatic index evaluates variation of principal stresses (pressures) with 
respect to density and represents rigidity of the gravitating fluid. We consider 
F to be constant throughout the stability analysis of the fluid. Equations (l22|l 
and (1261) lead to 


Pr 

P<t> 


-F 

-F 


b Cpo+ P<po , Pro ^PO , ProAo(j)o - ‘ 

-^PrO + -^- Pro + -^-— + -^ 

-Do ^ Po + PrO Po + PrO 00 PO + PrO 

b p0+ Pro , Cp^O , P^O^’PO , P<AO^O0O 

-Do Po + PfpO [po + P(/)O)0O Po + P 0 O 


T, (27) 
T.(28) 


Using Eqs. flSTj) . fl22|) . fl27|) and fl28|l in fl23|) . we construct a hydrostatic equa¬ 
tion given by 


Pro 

(Po 


^ (Po + P<t>o) cT 1 
Bo (Po + Pro) r Po + Pr 


-Ao(j)oKiT 


P<t>o) cT 


[po 

r 


Pro) 

1 

Po T P(j>o 


xAo(j)oK,T]] 


A' 


PrO 

Po + Pro 
Ao(poKiT 
bT 


Aq(I)qKiT 


Pr 


+ Mr 


P<t>o 


FT TfeT 

Pro yw 

L L-^o 

bT [po + Pro) 
Bq (po + p<^o) 
1 


^ , (Po + P0o) cF 
5, ' ■ ' 


A' 

x/-F p, 

^0 

+0O-Bo-^2 = 0. 


^0 

bT 


[PO + Pro) Po + Pro 

cT 


D (PO + Pro) + (PO + P 0 o)- 1 “ AQ(j)QKiT 

Bo r 


^ (Po +P0o) cT ^ 


Bo (PO+Pro) PO+Pr 


-Ao(j)oKiT 


00-So 


(29) 


This represents the general form of collapse equation which describes the 
instability of hydrostatic equilibrium of gravitating hlaments in BD gravity. 


3.1 Isotropic Fluid 

Here, we analyze instability ranges of isotropic fluid in Newtonian and pN 
limits. In isotropic fluid, all principal stresses are equal [pr = Pt/, = Pz)- Using 













































this condition in Eq. fl29p we obtain the corresponding collapse equation 

r 


" 

\bT cT 1 , , - 1 

- f 

" 

Pro 

n + 1 Ho0o-^i-^ 

-T 

Pro 


[Bo r Po+Pro \ 




bT cT 


+ 


PO + Pr 


-Ao<poKiT 

+Ao(PoKiT] 4^ - r 
+ (PoBqK2 = 0. 


^0 


bT cT 

{PO+Pro){-^ + —) 
-Do ^ 


Pr 


bT cT 

TW + - + 


-B, 


0 


r Po+Pr 


-AqcPqKiT 


4>oBo 


(30) 


Newtonian Limit 

The Newtonian limit in BD theory leads to the following 


mo 


Po s> Pro, Po s> p</, 0 , -Bo — 1, ^0 — 1 — 

(j)o = constant, Vq = V = 0. (31) 

Using these limits along with fl25|) . the collapse condition turns out to be 

r 


{PtoZn),! 


r'^c^ 


Po^N—^ + -^3 < 0, 


r 2(.2 


which gives 


r < 


PqZn 


mo 


-K^ 


\{proZ n) ,r ~jS^ProZN\ 


(32) 


Here 


Zn — [b + 


K. = - 


a'mo $ mo $ 

+ y—(1 A — ^PrQ — 


r‘^c^(j)Q (po 


This shows that the adiabatic index depends on dynamical properties such as 
density, pressure, scalar field. To preserve difference between configurations 
of pressure gradient and gravitational forces, we assume T > 0. Thus the 
celestial objects remain unstable until (1521) is satisfied which leads to 


PoZn^ - Ks 


\_(^ProZ n') ,r j. 2 ^ 2 ProZN~\ 

This leads to the following possibilities: 


> T > 0. 


(33) 
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1 . — K 3 = [{proZN),r —^^ProZ]\j\] 

2. PoZj\f^^ — < [(prO-^Ar),r —^^ProZjq\] 

3. PqZ]^-^^ — > [(prO-^Af);r 


The first and second case along with (j3^ show that the isotropic system 
becomes unstable for 0 < T <1. The corresponding expressions lead to 



(34) 


(35) 


These are the constraint expressions for a collapsing cylindrical isotropic 
hlamentary structure with 0 < T < 1. In the third case, the denominator 
is less than its numerator and hence in fl33|) . T can be taken greater than 1. 
The corresponding instability constraint is given by 



(36) 


for which T > 1 and isotropic cylindrical system becomes unstable. It is 
obvious that if the system is unstable for T > 1, then it will also be unstable 
for 0 < T <1. 

Post-Newtonian Limit 

The pN regimes are found upto order c~^ by taking 



where 


a = 


1 + ^BD 

2 + oJbd 


Using pN limits along with Eq. fl2^ in fl2^ . we obtain 


[{po+Pro)X,N]{^-2^)-K, 
[[ProXpN^ - [ProX,N] (^ - 2 ^)] ’ 


(38) 
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where 


XpN — 


(i(i _ ^!!!^) + £) + _^ 

rc^ r Po+Pro 


K 4 and are given in Appendix A. The expression 0381) describes condition 
for instability of a cylindrical filamentary structure in pN limits. Similar to 
the Newtonian case, the system collapses for 0 < T < 1 with the following 
constraints 


1 2r(^-^)+r Y„Mdr' rr _ T Y„Mdr' 

1 . PrQ = Xpj^e ■'’■0 


'■pN 

ipoi^-"^) + X, 


pN 


</>o 


(1 + 7^(1 - a)) + ^(1 + 4a) ) - 


( mo _ 4m2 2aml ^ 

^ 3^4 “r „ 4^4 ) ? 






(Po(5^ - Sj ) + yy 1^(1 + ^(1 - a)) + ^(1 + to )) - yy 7 |(., 


mo _ 4mg 2amg ' 

/p2^2 y»3^4 I ,p4^4 y 


—A ,2 


In the third case, T > 1 leads to unstable configuration with the following 
constraint 




'■pN 

mo 2ml , 


'ro 

2 mi 


TUn 


(1 + ~ “)) + 


rc^ 




mo 


4mn 2amn 




+ 


^ 4^4 


), 


where Yp^ = X. 


-1 

'-pTV 


^(1 + ^(1 - a)) + 73^(1 + 4a) . In this case, 0 < T < 


1 is also an instability range. 


3.2 Anisotropic Fluid 

Here, we have pro 7 ^ p^po 7 ^ Pz and hydrostatic equilibrium is described by 
Eq.([29]). 
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Newtonian limit 


Using Eq. fIST]) in ([29]), we obtain condition for unstable anisotropic filaments 
as 


0 < r < 


2pr0^ (PrQ—P^o) f cl 

(po ' r L r J 


+ Po{Zj\f + Kq) + 


IPrO^A-l' + - PrO^A-S? 


(39) 


where Kq and Kj are mentioned in Appendix A. Similar to isotropic case, 
this implies that for 


.J rO 


Pro < 

+ Po(l + A’e) + Ky] dr'] , 
r lies in (0,1) and the system collapses. If 


rZ7ve“^^“-^’'o l+Mr] 


Pro > 


( rT-o r *1 1 r£iu,,''i 


.J rO 


P4>o 
L r 


+ Po(l + ^o) + K 7 ] dr'] , 
the system becomes unstable for F > 1. 


Post-Newtonian Limit 


The collapse condition of anisotropic cylindrical hlaments in pN regime is 


r < 


(Pro + Po) [UpN] - K, 


[ [PrO 


pN\ 


+ ^-fu, 


pN 




where 


UpN 

VpN 


W. , b0O + Po) c , 1 

“rc2 {Pro + Po)r {Pro + Po) j’ 

rc^ (P 0 O + PO) r (Pfl)o + Po) J 


(40) 


The values of and Kg are given in Appendix A. The system becomes to 
unstable for instability ranges 0 < T < 1 if 
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• PrOrUM^-^)-{PrO+Po) [Upn]-K^ = [[ProUpN^ + ^-fU^N " ^V^piv] 

• PrOrUj,N{^.-^)-{Pr,+Po) [^piv]-i^9 < [\ProUpN^ + " ^^piv] 

and becomes unstable for F > 1 if 

• Pr0rUpN{^-^)-{Pr0+P0) [UpN]-K^ > [[ProUpN^ + - ^f"piv] 

4 Concluding Remarks 

The study of structure formation in modified gravity is an important issue. 
Cosmic filamentary structures with cylindrical symmetry arise in astrophysics 
both on large (cosmic web) as well as small (tidal tails) scales. The behavior 
of these structures has an important role in the formation of structures in the 
universe. In this paper, we have investigated instability ranges of anisotropic 
cylindrically symmetric collapsing hlaments in BD theory. We have used 
contracted Bianchi identities to obtain two dynamical equations of collaps¬ 
ing hlamentary system. By applying perturbation technique on BD as well 
as dynamical equations, we separate the unperturbed (static) and perturbed 
(non-static) distributions of all dynamical relations. We have developed hy¬ 
drostatic equation (collapse equation) through perturbed conhguration of 
second dynamical equation. 

The equation of state involving adiabatic index controls the ranges of 
instability for a collapsing hlamentary structure. We have used collapse 
equation along with equation of state to investigate the instability ranges of 
both isotropic as well as anisotropic BD huid at Newtonian and pN limits. It 
is concluded that in both approximations the adiabatic index depending upon 
dynamical properties (energy density, pressure, scalar held terms and some 
constraints) controls the instability ranges. We have constructed constraints 
on static radial matter pressure under the ehects of scalar held. It is found 
that the cylindrical hlamentary structures always remain unstable for 0 < 

T < 1, while T > 1 is the instability range for the special case. We would like 
to mention here that the instability ranges for spherical as well as cylindrical 
distributions in GR depend upon T < I and T < 1. In f{R) and f{T) 
theories, physical variables such as density, pressure and respective modihed 
dark terms provide the instability ranges. The instability range of spherically 
symmetric isotropic BD huid is r > I while anisotropic spherical BD huid 
always remains unstable for 0 < T < 1 and T > 1 leads to collapse only for 
the special case. 
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Appendix A 

The non-zero components of BD equations for interior region are 


G,o = \iT^o+Ti) = \(pA^ + ‘^ 


2 , 


yl20'2 

52 


2A B (P'A^ /B' C'\ A^cj)" 

(P\lA^B^cj^ \B^~c)^ B^ 20 ^ 

^ 1 ^ ‘^BD 50' 

Goi = -(5oi + Toh) = -^(00 ) + T 


A 


5 


(41) 

(42) 


2 J,2 


Gn = + rS) = ^ I p,B^ + + !^) I + 


B‘^ct) I A C 

[a^c 


20 

b' (A C'\ B^V{(j)) 

'A^~c)^ 20 ^ 


A^ 


B(f) 

A^ 


UJBDpC^^^ C^'\ 


(43) 


G 22 = -Xt^2+t!,) = ^ I + ^(^ - ^) 1 + 


0(^2 


20 ' 212 52 M ^20 

(A B\ C^(l)’ (A B'\ C‘^(i)" C^V{(j)) 


+ 


+ 


2120 \ A Qj B^(f) \A Bj 520 

a2 j/2 




A 2 52 


+ 


+ 


20 


(44) 


2120 

and Eq.(l3]) becomes 


ABC 
I ^ 5 ^ C 


520 


A' 5' C" 


2420 

"" (45) 


520 20 


24 


5 




2 l A25 2125 

1 




A' 5' C" \ 0" 

~C^) ~ 


‘^ojbd + 3 




(p -2 3pr +p<j>+ Pz) + ~ 24/ 


( 46 ) 


The scalar terms Ki and K 2 of Eqs. flT^ and flT^ are 


Ki = Tot, 21 


t-i 


-01 


A-^B-^ + 


-1 


+ Tot,2l 


-2 \ i-2. 
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+ - 2T^^B^AB' - T^^A-^B-^A', 

K 2 = Tf^B-^(t)'(j)-‘^ + {^p + T^^ A-‘^B-^(t)-‘^^- {t^^A-‘^B-‘^'^ B 
- {T,,B-^),^B. 


The static distribution of BD field equations is 


2B14-1 


Po , i^bd4>^ 
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Pro ^BD(pO 
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_^0 ^0 Bq 

00^0 00-^0 
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AqBIcPo BlrcpQ 

^000 

^0-^0 00 
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1 A?, 


200 Blr A 

00 ^0 

5q0o 200 




00 -^0^000 


00 ^BDfp', 


^0 00 
1 


i) 

Bi 


/2 


dloi^o" 


Blrcpo 
S3 U Ao 


+ 




200^0 200 


(47) 

(48) 


(49) 


(50) 


The unperturbed wave equation is 


0Ao 00-^0 




R2 

-°o 


+ 


r So 2 ujbd + 3 


[(PO + 3pr0 + T(/)0 + Pzo) + (00^0 — 2Vo)] . 


The static distribution of Eq. flT^ is identically satisfied in static background 
while (IT^ turns out to be 


Pro 

-Bo 00 
where 


, 0oP»’O , r , N ^0' 
+ ,o + (Po + Pro) 


^iBi 


+ 


AoSo0o So0or 


■(PrO P0o) 


A = 0, (51) 


K' — T 

^2 — llunp 


4 00 _ PrO I ^ A (p^o PT I -Bq (Bii^^p); r 

''“""0oB2 0oSo ^^^°+^°AoSo0o Bor0o + So 


and the term Tfi«np represents unperturbed form of energy tensor due to 
scalar held. 


15 


































The static part of Eq. ffTII]) is 


s(=) 

PrO 
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(52) 


The perturbed form of BD held equations are 
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and perturbed wave equation is given by 
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Here T^t(nnp) and indicate unperturbed as well as perturbed distribu¬ 

tions of BD energy part, respectively. 

The values of K 4 and in (155]) are given by 
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The values of Kg and Kg in flT0|) are described as 
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